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The Galactic Cepheid period histogram has a strong dip between 8 and 10 days that has defied an explanation based on 
evolutionary and linear pulsation studies. We show here that this deficiency is caused by the instability of the nonlinear 
fiindamental pulsation cycle in this period range. The strong metallicity dependence of this instability is consistent with the 
absence of a corresponding minimum in the Magellanic Cloud data. Our results also suggest that the Galactic Cepheids must 
have a large spread in metallicity. 
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It is a well known fact that the observed period distri- 
bution of the Galactic Cepheids has a pronounced mini- 
mum in the 8-10 d period range, as shown in Fig. 1 which 
displays the histogram for the period distribution of the 
Galactic Cepheids constructed from the Galactic Cepheid 
Database of Fernie et al. (1995). Since it is difficult 
to separate overtone and fundamental pulsators for the 
Galaxy above 5 d our histogram necessarily contains both 
the fundamental and the overtone Cepheids. 

Becker, Iben & Tuggle (1977) combined the location 
and duration of Cepheid model crossings of the instabil- 
ity strip from evolutionary calculations with a birth-rate 
function to infer a theoretical period histogram. One 
of their conclusions was that a minimum in the distri- 
bution was not compatible with a standard birth-rate 
function, and could only be explained if an ad hoc two- 
component birth-rate function were adopted (cf. how- 
ever Chiosi 1989). Their evolutionary computations were 
performed with the now superseded Los Alamos opaci- 
ties which are now known to be considerably too weak. 
However, while the new opacities will produce a different 
period distribution it is difficult to see how they might 
cause a two-humped one. 

In this Letter we show that the minimum in the period 
distribution of Cepheids is a result of the nonlinear dy- 
namics associated with the fundamental pulsations of the 
Cepheids, and that it has therefore nothing to do with 
evolutionary calculations. 

The proper mass-luminosity (ML) relations for the 
Cepheid model sequences would normally have to be 
obtained from evolution calculations. However, at the 
present time there is enough uncertainty and disagree- 
ment (Buchler et al. 1996, Beaulieu et al. 1997) to lead 
us to determine these ML relations differently. From the 
structure of the Fourier decomposition parameters 02 1 
and i?2i for the Galaxy and for the Magellanic Clouds 



one infers that the 2:1 resonance occurs in the vicin- 
ity of a period of 10 d. Based on this fact we con- 
struct the ML relation as follows: For each metallic- 
ity Z we determine the mass M and luminosity L for 
which the Po=10d equilibrium-model is resonant viz. 
^2/^0=1/2 (bump Cepheid) with an effective temper- 
ature Teff that lies AT^IOOK degrees to the right of 
the fundamental blue edge. (Mathematically, for given 
composition parameters X and Z and a given AT we 
solve for Po{M, L,Teff)^10, P2{M, L,Teff)=5 with Teff 
= Tbe-AT where Tbe ^ Tbe{M,L)). From these an- 
chor values (M, L) we then derive a ML relation with a 
slope chosen to be 3.56 that is close to what evolutionary 
calculations indicate. 
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Fig.l: Period distribution of fundamental and overtone 
Galactic Cepheids d'apres Fernie et al. 
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We use the Livermore OPAL95 opacities (Iglesias & 
Rogers 1996) combined with the molecular opacities of 
Alexander & Ferguson (1994), to compute the nonlin- 
ear fundamental pulsations with the relaxation method 
of Stellingwerf (1974). We put a temperature anchor 
at 11,000 K with 30 constant mass shells in the sur- 
face, and with 90 geometrically increasing zones inward. 
The pseudo- viscosity parameters are Cq=4 and a=0.01. 
Convection is ignored, and we caution that the models 
lose their validity far from the blue edge. Concomitantly 
with the relaxation to the periodic pulsation the code 
performs a Floquet analysis of the limit cycles (i.e. the 
periodic finite amplitude pulsations) (e.g. Buchler 1990). 
We recall that the Floquet exponents measure the linear 
stability of the limit cycle to perturbations with respect 
to all the possible modes (e.g. Ince 1944): For a limit 
cycle to be stable, and thus to be observable, all Floquet 
exponents, A^, have to be negative. This stability anal- 
ysis is an extremely useful byproduct of the relaxation 
method. In fact without the Floquet analysis we might 
have to integrate thousands of cycles to ascertain stabil- 
ity of a limit cycle, and, in the case of instability, we 
woiild not be able to compute the limit cycles at all and 
determine how unstable they are. The relaxation code on 
the other hand is robust enough to converge on a limit 
cycle even when the latter is mildly unstable. 
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Fig.2: Top, middle: Floquet stability exponent Ai for 
the nonlinear fundamental Cepheid pulsations as a func- 
tion of period and metallicity Z (with X=0.70) for mod- 
els 100 K to the right of the fundamental blue edge; Bot- 
tom: Square of the radius Fourier amplitude of the limit 
cycle scaled by the period, ap, cf text. 



In Fig. 2 we exhibit the behavior of the Floquet sta- 
bility exponent, Ai, of the fundamental limit cycle as a 
function of its (nonlinear) period PqK (In this period 
region it is the perturbation with the first overtone that 
is the least stable). The six sequences have metallicities, 
with Z ranging from 0.013 to 0.035. In order to avoid 
cluttering we have split the figure into two subfigures, the 
top showing the metalhcities Z= 0.013, 0.015 and 0.020, 
and the bottom the values 0.020, 0.025, 0.030 and 0.035. 
The hydrogen mass fraction has been chosen at X=0.70. 
All model sequences run AT=100 K to the right of the 
fundamental blue edge. For metallicities in the range Z 
PS 0.014 to 0.035 these fimdamental limit cycles are thus 
unstable to a perturbation in the first overtone. 

In order to see the effect of the location of the Ccpheids 
with respect to the blue edge we display in Fig. 3 results 
obtained for a sequence with AT = 400 K (with the same 
ML relation as for the 100 K sequences described above). 
The width of the period region of unstable limit cycles 
thus depends somewhat on location with respect to the 
instability strip, and shrinks with AT. 
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Fig. 3: Floquet stability exponent Ai for the nonlinear 
fundamental Cepheid pulsations as a function of period 
for two sequences with AT =100 and 400 K and with 
Z=0.020 and X=0.70. 

Older nonlinear Cepheid pulsations that were per- 
formed with the now obsolete Los Alamos opacities 
(Moskalik & Buchler 1991) indicated that the fundamen- 
tal Cepheid pulsations were stable in the vicinity of 10 
d (except for an absolutely minute range, as their Fig. 1 
shows). It is quite clearly the overall increase in the opac- 
ities that is responsible for the instability. One can get a 
good feeling for the sensitivity of the Floquet exponents 
to opacity from Fig. 3 of Buchler (1996) where the results 
of calculations with the OPAL95, the OPAL93 and the 
Los Alamos opacities are compared. Roughly speaking, 
the change from OPAL93 to OPAL95 is equivalent to an 
increase of Z from 0.02 to 0.03, for example (Buchler, 
Kollath, Beaulieu & Goupil 1996). 

We caution the reader though that there is some sensi- 
tivity of the Floquet exponents to the zoning and to other 
numerical parameters that are used in the nonlinear cal- 
culations (Kovacs 1990, Yecko, Kollath & Buchler 1997), 
as well as to convection whose effects we have ignored. 
The precise values of the metallicity Z for the onset of in- 
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stability and of the period range thereof should therefore 
not be taken too literally, although the existence of an in- 
stability in the broad 8 10 d range seems to be essentially 
independent of such numerical parameters. 

What is the reason behind the instability of the fun- 
damental limit cycles? We show now that the resonance 
between the fundamental mode of oscillation and the sec- 
ond overtone {P2/PQ = 1/2) cause an overall decrease of 
the pulsation amplitude that in turn destabilizes the fun- 
damental limit cycle. We note that this is the same res- 
onance that causes the familiar Hertzsprung bump pro- 
gression. 

It is a well known observational fact that the pulsation 
amplitudes of the classical Cepheids exhibit a substantial 
drop for periods in the vicinity of 10 d. The same feature 
is found in numerical modelling. The amplitude equation 
formalism (e.g. Buchler 1993) explains how nonlinear ef- 
fects associated with the resonance P2/ Po=l/2 between 
the fundamental and the second overtone are responsible 
for the drop in the overall pulsation amplitude. Further- 
more, large amplitudes increase the stability of a limit 
cycle pulsation, as the same formalism shows. Indeed 
the Floquet exponent Ai behaves as (Buchler, Moskalik 
& Kovacs 1991) 



Ai = (ki - qwAl - 912^2) -Po 



(1) 



where Ki is the linear growth-rate of mode 1, q^j are non- 
linear coupling coefficients that depend on the structure 
of the star, and Aj are the pulsation amplitudes of the 
excited modes, i.e. the fundamental and the resonant 
second overtone here. 

For simplicity, let us make a few approximations. First 
we ignore the smaller amplitude of the second overtone, 
A2 compared to A^. Second, with good reasons we as- 
sume that the cubic coupling coefficients q are positive 
and that they scale as qio = ciqPq (Kovacs & Buch- 
ler 1989). Third we disregard the small variation of kiPq 
over the plotted range of periods. In Eq. 1 the amplitudes 
are relative, i.e. they refer to 5R{t)/R. Thus we associate 
Aq with the lowest Fourier amplitude oi of the absolute 
computed radius variations, scaled by the stellar radius 
i?*, i.e. with ap = ai/i?*. How the radius scales with 
period can be found by the following rough estimate. Us- 
ing the mass luminosity relation L ~ M^-^^, the period- 
mean density relation Pq ~ p"^/^ ~ i?^/^M~^/^, the 
surface luminosity equation L R'^T'^ and the shape of 
the blue-edge T ~ L", with a ^ —0.05, one finds 



^^^^^ ' ^-^"3.56(1^- 4a) 
With /3 « 1 Eq. 1 is reduced to 

Ai « Ki Po - cio ap 



1.27 



(2) 



A comparison of the bottom panel of Fig. 2 with the up- 
per panels clearly confirms the correlation between pul- 
sation amplitude and stability. 



On the basis of the stability properties of the funda- 
mental limit cycles exhibited in Figs. 2 and 3 we thus 
reach the following conclusions. First, for mctallicities 
in the range Z « 0.014 to 0.035 the fundamental limit 
cycle is unstable to a perturbation in the first overtone. 
Second, the window of instability shifts to higher period 
with increasing Z. 

Turning now to the astronomical implications we note 
that the Galactic period histogram has a dip, but not 
an actual gap in the 8-10 d range as it should have if 
the metallicity of all Cepheids were as large Z=0.02. 
If the minimum of the histogram is now interpreted as 
a dispersion in metallicity the observed w 8-10 d fun- 
damental Cepheids, approximately a third, must there- 
fore have lower metallicity, Z < 0.013 (or else an un- 
likely Z >0.035). We note though that Becker et al. 
(1977), albeit on totally different grounds, also suggested 
that Galactic metallicity dispersion is large (a range 

Zmax /Zfnin of 3 tO 5). 

However, it may be objected that since the stability 
of the fundamental limit cycles depends on the location 
with respect to the instability strip a wide strip could 
also reduce the Cepheid deficiency especially on the lower 
period side. There are evolutionary arguments (Becker 
et al. 1977) that the Cepheid instability strip may be 
much narrower than generally assumed. An estimation 
of the instability strip from the Fourier decomposition 
parameters also suggested a narrower rather than larger 
instability strip (Buchler, Moskalik & Kovacs 1990). 

How much of the survival of 8-10 d Cepheids is in- 
deed due to a dispersion in metallicity rather than due 
to the width of the instability strip can be tested with 
an observational measurement of the metallicities of the 
individual Cepheids. 

While there is a deficiency of fundamental Cepheid pul- 
sators in the 8-10 d range, there is absolutely no reason 
to believe that there is a deficiency of the corresponding 
stars. As expected, hydrodynamical calculations show 
that these stars pulsate in stable first overtone limit cy- 
cles. They should therefore show up as an excess in the 
histogram at lower periods (lower by a factor « 0.7, the 
typical period ratio Pi/Pq). Indeed, the histogram of 
Fig. 1 is certainly compatible with an excess in the 5.6- 
7.0 d range. (The left overbar shows the first overtone 
period range corresponding to the fundamental range 
shown on the right.) 

Turning now to other galaxies, we recall that Becker 
et al. (1977) present also a period histogram for M31, a 
galaxy that also has a relatively high metallicity. Again 
the histogram shows the dip in the 8-10 d range, consis- 
tently with our results. 

In contrast, the Magellanic Clouds are known to have 
considerably lower average metallicities than the Galaxy, 
and indeed, in agreement with our stability analysis, the 
period histograms for the Magellanic Clouds (Becker et 
al 1977) do not show much indication of a minimum. 
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The small hollow near 9d, if statistically significant, may 
again be an indication of a spread in metallicity putting 
some Magellanic Cepheids above the threshold value of 
Z « 0.013. 

If the 2:1 resonance at 10 d plays such an important 
role one may wonder if other resonances might make 
themselves felt, especially because the Galactic Cephcid 
period distribution (Fig. 1) seems to indicate a deficiency 
of Cepheids in two other places. 

First, if the gap near P « 25d is indeed significant 
it has an interesting implication. Numerical hydrody- 
namic calculations (Moskalik & Buchler 1991 and Moska- 
lik, Buchler & Marom 1992 for the new opacities) show 
that the Po/Pi = 3/2 resonance can destabilize the fun- 
damental Cepheid pulsation and lead to a periodic pul- 
sation with double period, in which alternating cycles 
differ only slightly however. We note though that Fer- 
nie (preprint) has not found any evidence of alternations 
in his Cepheid data. However, in support of this sce- 
nario Antonello & Morelli (1996) has suggested that the 
star CC Lyr does exhibit alternations. The observational 
data of CC Lyr are very limited though, and additional 
observations of this star are necessary to confirm the ex- 
istence of alternations. Is it possible that the dip in Fig. 1 
could arise because such stars with alternating cycles may 
not have been classified as Cepheids? If this resonance 
is indeed strong enough to cause instability to alterna- 
tions then this would add an additional observational 
constraint for Cephcid models. 

Second, Fig. 1 perhaps also suggests a minimum in the 
vicinity of 19 d. One notes that this is the period re- 
gion where the fundamental mode is in a 4:1 resonance 
with the fifth overtone {P5 / Po=l / 4:) . However, the linear 
stability of the 4th overtone is quite large, and a hydro- 
dynamical survey is necessary to verify whether it can 
play a sufficiently large dynamical role to destabilize the 
fundamental limit cycle. 

In conclusion, the simple and natural explanation for 

the deficiency of w 8 10 d Cepheid variables is that the 
fundamental mode limit cycles are unstable. The corre- 
sponding stars pulsate in the first overtone with period 
Pi « 0.7Po; giving rise to a relative excess of overtone 
Cepheids in the corresponding range of w 5.6 7.0 d peri- 
ods. It is thus not necessary to invoke an ad hoc two- 
component birth-rate function to explain the Cepheid pe- 
riod distribution. 

Our calculations and the agreement they provide with 
the minimum in the observed Galactic Cepheid distribu- 
tion, and the lack of a minimum in the lower metallicity 
Large Magellanic Cloud data also provide a further con- 
firmation that the new opacities are in the right ball-park. 
Indeed, if they were much weaker, the 8-10 d fundamen- 
tal Cepheid pulsations would all be stable, thus not giv- 
ing the observed minimum. On the other hand if they 
were much stronger, then they would predict a minimum 
for the LMC as well, and a smaller one, or none at all for 



the Galactic Cepheids, in contrast with observation. 
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